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Abstract 

We show that formally for the standard ansatz relativistic point-coupling mean-field 
(RMF-PC) model a lagrangian density C is not equivalent in Hartree and Hartree- 
Fock approximations. The equivalency can be achieved only if we use a "complete" 
ansatz at the cost of introducing new parameters in the model. An approximate 
treatment of the exchange terms from standard RMF-PC indicates that these effects 
cannot be easily, if at all, absorbed by a Dirac-Hartree approximation. 
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1 Introduction 



Relativistic mean-field models, which describe the nucleus as a system of Dirac 
nucleons that interact with each other by exchanging mean meson fields, have 
been successful in describing nuclear matter and ground state properties of 
finite nuclei [1,2,3,4,5,6,7,8]. Applications include the calculation of deformed 
nuclei [8,10,11,12], odd nuclei [8], nuclei at the driplines [13], the potential en- 
ergy surfaces of heavy nuclei [8], and the prediction of superheavies [9,8,14]. 
Another variant is the point-coupling model [18,15,17]. The difference be- 
tween this model and the Walecka model is the replacement of the mesonic 
potentials of the Walecka model by explicitly density- dependent potentials. 
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Nikolaus, Hoch, and Madland [17] used the Hartree form of this model to 
calculate some observables of finite nuclei and nuclear matter. They obtained 
similar predictive power as the Walecka model, but used different weights and 
observables from Ref. [1] to obtain their parameter set. In Ref. [18], Rusnak 
and Furnstahl have shown the profitability of applying the concepts of effective 
field theory to the point-coupling model. A more recent discussion about the 
point-coupling model can be found in Ref. [15,16]. It consists of a careful way 
to obtain a new parameter set for the point-coupling model which is likely to 
correspond to the optimal minimum by combining different methods for the 
X 2 minimization. This parameter set was applied to a wide area in nuclear 
theory, from spherical finite nuclei, axially deformed nuclei, superheavy ele- 
ments, potential energy surfaces, nuclear matter, neutron matter, up to exotic 
nuclei. It could be shown that point-coupling models can deliver a predictive 
power for finite nuclei comparable to models employing meson exchange. On 
the other hand, especially concerning isovector properties, the point-coupling 
model suffers from the same unresolved systematic deviations. 

From the success of the RMF (Walecka and point-coupling) models, we sup- 
pose that the exchange effect is already absorbed effectively in the coupling 
constants of the model due to the fitting procedure, but some calculations 
which take into account the exchange explicitly by using the linear Walecka 
model have shown that this does not seem to be the case [19,20]. Understand- 
ing this situation seems necessary and interesting. Unfortunately, due to the 
finite range of the mesons and the mesonic nonlinearity of the RMF model, 
it is difficult to expand this model in the above direction. This problem does 
not appear if we use the point-coupling model [17], because in this model the 
linear and nonlinear terms are explicitly density dependent. Therefore in this 
paper, we choose the point-coupling model for a first study of the role of the 
exchange effects. 

We would like to stress that we regard this work as a model study. The aim 
is to check the principle effects of exchange terms as compared to an effective 
Hartree theory. This is done for the two-body point-coupling terms of the 
standard RMF model. Higher order couplings are treated without exchange 
and the pseudoscalar (pionic) channel is omitted for the time being. 

In the following, we will use the phrase strict Hartree for an ansatz that is 
taken from the start in the Hartree approximation. We will use Hartree-Fock 
for a model in which the exchange terms are treated explicitly. If this model 
can be transformed into a form consisting of local currents and densities only 
(no matter how complicated they are), we will call the result of this rewriting 
an effective Hartree theory. Mapping formally this general form onto a Dirac- 
like structure leads to the Dirac-Hartree theory. 
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2 The formulation of Hartree-Fock Point-Coupling Model 



2.1 Basic one-particle operators and densities 



We beginn by defining the basic Dirac & isovector matrices F,-. There are two 

sets 




H ={l,7^,7 M r} 

HF = H U {f , cv, cr^f , 75, 7 5 f, 757^, 757^} 



The "Hartree" set H covers the currents initially given in the RMF-PC La- 
grangian. The "Hartree-Fock" set HF covers all conceivable Dirac matrices as 
they are produced by the Fierz transformation [28]. The currents are defined 
as: 

where 8^=8^ — <9 M . These are the current operators from which actually only 
Jj is used at operator level. The associated currents are the expectation values 



7" . — r .\i/ A - — 8 7 

a 

<— > 

**=E*«tr r ;*« . c, = £wy^ a , (i) 

a Z< a 

where the sum is running over occupied particle states only (anti-particles are 
neglected = no-sea approximation). 



2.2 The Lagrangian 



Starting point is the Lagrangian 



^= 9 EM 2 + H 2 ]» (2b) 
z jeR 

^1=\psJl + \lsJ A s . (2c) 

At this point, replacing the current operators by their expectation values leads 
to a strict Hartree theory, as is mostly applied to nuclear structure problems. 
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But we are going here for a Hartree-Fock treatment. To this end we chose 
a Slater determinant |0) as the ground state with A occupied single-nucleon 
levels (in the positive energy sector). It is defined by 



V i < F, 



(3) 



(F indicates the fermi surface). We compute the expectation value of the La- 
grangian (2) in standard manner and employ Fierz transformations to express 
exchange terms through simple currents. This yields 



(0\£?°\0)=££ C -£Z 



PC 



-PC 



rPC 



=E 

ieH 



\ ^ J? + \ k Aj 



r PC 



EE 

igH j'eHF 



C 



B] + J, V 



(4a) 
(4b) 

(4c) 



The constants Cji stem from the Fierz transformation, their explicit form can 
be found in table 1. 



Thus we have performed - by virtue of the Fierz transformation - a map- 
ping from a Hartree Fock model to an effective Hartree model for all linear 
terms. The expectation value of non-linear term has been calculated in 
Refs. [27,15]. The results are quite lengthy, therefore it is not useful to consider 
them explicitly here. The exchange terms emerging from the nonlinear part 
of the Lagrangian form a large variety of density-mixing terms, involving also 
products of isoscalar and isovector densities. This leads to modifications in 
both the isoscalar and the isovector channel of the effective interaction with- 
out introducing new parameters. Another feature of this result is that the 
contributions from the nonlinear exchange terms cannot be absorbed fully by 
a redefinition of the coupling constants into terms of an d they also cannot 
be considered as small corrections. 



2.3 The Hartree-Fock equations 



Variation of the effective Lagrangian density (4) with respect to the single- 
nucleon wavefunctions yields the effective mean-field equations 
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Table 1 

The Cji constants of equation (4). 






= -il^ + m B + V j Y 3 + U^Tjdn + W^jd^] V a 


(5a) 


v% 


= s + la v a 


(5b) 






(5c) 




= W l + -f a W 2a 


(5d) 


s 


= a a J a - (5 a + ~5 a )d^A^ a + p a (id^ a + C a ) 


(5e) 


ya 


= a b J b a - (8 b + ~8 b )d»A« b + + 


(5f) 




4 4 


(5g) 


w 1 


2 ~ 2 ~ 

= — -S a J a , Wi a = — -SbJba- 


(5h) 



A double index a/b implies summations over S,D/V,R. The constants in front 
of the densities in Eq. (5h) (a a / b , 5 a / b , 5 a /b) are functions of the meson coupling 
constants and masses. The 5 a / b , S a / b constants reflect the combinations of Fierz 
constants Cy, the meson coupling constants and masses of the mesons from 
the derivative exchange terms of Eq. (4) in the equation of motion. Because 
the formulas for the constants are lengthy and do not give particular insight, 
we do not show them. S, V a , W 1 , W 2a are of order < 1 and U^, U afl are of order 
v . We can see here that the Cj densities give rise to the term with prefactor 
W j and the B ifl the term with prefactor in the equation of motion. 

Equation (5a) starts out like a Dirac equation and continues with very uncon- 
ventional terms oc Yjd^ and oc Yjd^d^. The latter are generated from the Fierz 
transformed exchange terms oc B and C in Eq. (4). These terms are awkward 
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to handle. We would like to transform Eq. (5a) into a more standard form. 
This will be done in section 4 by means of a Gordon decomposition. Before 
that we first check the importance of the new terms. This will be done in the 
following section. 

From the discussion above it is clear that, in contrast to the claim of Ref. [17], 
qualitatively, if a standard ansatz is used, a Lagrangian densitiy C HF can 
be determined in a relativistic Hartree-Fock sense that is not equivalent to 
that determined in a relativistic Hartree sense, C H . Of course, the Hartree or 
Hartree-Fock calculations belong to effective theories, where many details or 
ignored corrections can be hidden in the form of the effective interaction and 
in its fitted parameters. In this sense, the models can only be compared after 
refitting the parameters. Therefore, only quantitative calculations can actually 
distinguish the predictive capabilities of the models. Before embarking into 
the actual calculation, though, we present a rough approximation that gives 
physical insight into the exchange contributions. 



3 Estimating the Densities £> JM and Cj 

For a first overview, we compare the order of magnitude of the energy density 
contributions of each term in Eq. (4) in the plane wave (nuclear matter) limit. 
The dominant contributions in the energy density come from the terms with 
Fj — 1, 7 . The terms with Tj = 7 show a clear order counting and are easy 
to calculate. It is illustrative to examine their order of magnitude. They will 
be represented by : 



ei = 


--a v J v 


£2 = 






-—a v J v 


e ex2 = 


c v 5 v B v 



»-T^£[*«7o(Vtt Q ) - (W a ) 7 o^J.[^/37o(W /3 ) - (V^) 7 o^] 

4 a/3 
€ lx2 = c v$vJv Cy « C v Sy 

a/3 

c v is a constant due to the Fierz transformation. It is a function of the coupling 
constants and fulfills c v < 1. Because the formula for c v is lengthy and does not 
give particular insight, we do not show it here. In the last two equations {e\ x2 
and e 2 ex2 ), we neglect the energy dependent parts (retardation parts) by as- 
suming that this contribution is small. Then in the plane-wave approximation 



6 



Terms 


pynrPKSi on 




order 


c i 


1 ~2 


Hirprt, nonrlprivativp t.PTm 


v° 


f o 


~ 


dirpct dprivativp tprrn 







~ -\cyavpQ 


exchange nonderivative term 


v° 


A 

t ex2 


r S 4 ~9 2 

Mvg^PoPf 


exchange derivative term 


v 2 


t ex2 




exchange derivative term 


v 2 



Table 2 

The order counting 

the order counting which is shown in table 2 is obtained. We use 

(2S+1)(2/+1) r j3 - 
Po = ^3 J d 3 kV kl0 ^ k 

where S and I denotes spin and isospin, respectively, and pf is the fermi 
momentum. If we define pp/Tn v = v, every exchange term in the derivative part 
gives a contribution of order v 2 , where v < 1. Then in the nuclear matter limit, 

Pi 7^ 0, V Pi =0, | Bi | /{rriBpi) is of order v and | Ci \ /{m 2 B pi) is of order v 2 , 
while in finite nuclei, we know that Pi ^ and Vp« t^O, where | Vpi \ /{rriipi) 
can be estimated to be of order v or maybe even smaller [17] by comparing the 
direct and derivative terms in the Hartree level of the point-coupling model. 
Here the index i denote V or S. If we consider the densities in finite nuclei 
as consisting of the nuclear matter value plus fluctuation corrections due to 
quantum effects, it seems reasonable to assume in finite nuclei that on the 
average | Bi \ /(rasPi) is also of order v and | Q \ /(jn 2 B pi) of order v 2 . 

In conclusion, it is reasonable to assume that the contribution from the deriva- 
tive exchange terms could be similar or even a bit larger than the direct 
derivative terms, so that they cannot be neglected if the direct derivative 
terms, which are necessary for a reasonable description of the nuclear surface, 
are kept. 

It should be noted that even the nonderivative exchange terms are larger 
than the derivative terms in order of magnitude, but their largest contribu- 
tions (scalar-isoscalar, vector-isoscalar and vector-isovector) can be absorbed 
into the nonderivative terms through a redefinition of the coupling constants 
and refitting of the experimental data. The crucial difference between Hartree 
and Hartree-Fock calculations of the nonderivative terms is represented by 
the scalar-isovector term, which is usually omitted in most Hartree models. 
In the Hartree-Fock approach presented here, it emerges naturally without 
introducing an additional parameter. 

As a further check of our procedure, we show a transformed form of the B a 
terms through the Gordon decomposition [29] by using the exact equation 
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of motion in Eq. (5). Here we only calculate the scalar case; the vector case 
is similar but not as simple. In the limit of small v 2 (for more details see 
Appendix (2)), we have 



■^m a ) - to*^ q7 ct # q - V^ a m 



a 



(6) 



a L 



The same result as above is obtained by using the equation of motion in the 
form of a Dirac equation with scalar potential (m*-m) and vector potential 
V a . It means that in the limit of small v the B a can be determined from the 
Dirac equation (second-order differential equation). 

A similar situation holds for all Ci and £>j M , so that we can say that the trans- 
formed form of the exact Ci and £>j M consist of parts responsible for generating 
effects beyond the Dirac equation. The order of magnitude analysis above 
shows that these parts are small. Neglecting them everywhere (via projecting 
onto a Dirac-Hartree structure) will lead to the Dirac equation but neglecting 
these small parts means also neglecting the retardation and nonlocal effects. 
This is what we will do in the next section, where we map the exact Lagrangian 
including exchange onto a Dirac-Hartree structure. 



4 The Approximate Lagrange Density 

The aim of this chapter is to derive an effective Lagrangian which leads to 
a Dirac equation without extra derivative terms, i.e., we will formulate, by 
employing some approximations, a Dirac-Hartree model. 

The approximate forms of the densities Ci and B^ are not yet known. We 
expect to obtain from them an approximate Lagrangian density (with the 
approximate forms of the densities Ci and B^), which should lead to a Dirac 
equation of the form 



where to*, V a and T a " are real functions. It implies a connection between the 
two types of functions: to*, V a , T a/3 and J7i, A^, & , Bj^. Therefore we will 
try to obtain the C j and B^ as functions of the m*,V a , and T Q/3 . This will be 
done in an iterative procedure starting from to* — to, V, T — and generating 
succesively terms with increasing orders. Finally, we use the approximate form 
of to*, V a , and T a/3 to calculate C J and B^ and from these approximate 
densities C J and B^ we obtain the approximate Lagrange density. These steps 
will now be done explicitly for the scalar case (r =1). 



+ i(m* + la V a + a af3 T al3 )]^ a = 0, 



(7) 
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Starting point is the Gordon decomposition [29] of Eq. (7) which leads to 



*ah^ ~ i(m* + l«V a + a aP T^)\ ffl! a 
V* Ah^ + i(m* + la V a + a aP T a ^ a = 0, 



(8) 



/i = a^^fj, and a M is an abritrary vector. Chosing fi — lp, yields 



B s = ^{^ a {d^ a )-{d^ a )m, 



a 



) 



Here f(T a p) denotes a complicated function of T af3 which is zero for T a/3 =0 
and the order of magnitude of this functions can be estimated as small there- 
fore it is not necessary to show this function explicitly These are the desired 
expressions for B$ and Cs in terms of m*, V, and T. We would like to express 
them in terms of standard currents. The construction for that proceeds as 
follows: 

The iteration procedure leads to m* = mB + asPs(x)+..., V u = ayJ v {x) +..., 
and T a/3 = where as and ay are effective coupling constants. They are 
functions of as, ay , 5s and 5y, which come from the contribution of the 
direct and the exchange part in front of the nonderivative and the deriva- 
tive of the scalar and vector densities, respectively. In the above equations, 
... denotes contributions from the derivative, isovector, nonlinear and tensor 
correction parts. They will generate nonlinear terms in the Lagrangian density 
that consist of all possible combinations of the above mentioned densities and 
of all possible order. Now we obtain the approximate forms by neglecting ... 
(neglecting higher order terms) 



B s = -i[d^ a a^ a ) + 2{m B + a sPs )^al v ^a + 2a v J v ^ a ^ a ], (11a) 
C s = [(m B + asPsY^a^a + 2(m B + a sPs )a v J v ^ all/ ^ a + l/2^9„(* a * a ) 



From Eq. (7) we deduce directly 



= [m* 2 * A + 2m*V* a >y v * a + l/2c>^(* A) 
+V v Vjl a ^ a - d v {V^ a a^ a ) + t a /(T a/3 )tt a ]. 



(10) 



yielding 



+ {a v J v ) 2 ^ a ^a - d v {a v J Vl ^ a a^ a )\. 



(lib) 



Putting the pieces together yields 
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-2m B J^{J T ^) - l -d»{J T ^)d a {J T ^), 



^s C £ )= ^Ba s p 3 s + 2a 2 s p% - Am B a s ps4 + 2a v p s J T v^ v {J») 



-2a sPs J^{J T ^) - 2a 2 s p 2 s J 2 v . (12) 



The densities and currents appearing in these formulas are special cases of the 
general Ji in Eq.(l), they read 



a a 
a a 

J7(z)=E^<r^a, *W = £* Q ^f* Q . ( 13 ) 



We see that £/ie approximate form of the interaction effects in C and corre- 
sponds to having all possible terms made up of products of different densities. 

Eq. (12) define the scalar case of the approximation for the derivative exchange 
terms in Eq. (4). The other cases can be derived similarly. The results are 
displayed in table 3. Inserting these results into C[, c cS = (0o| : £l° '■ I0o) 
yields: 



^Kf=E*«( < 7^-'»B)*« 

a=A 

~\ & sfis - \ & tspls - \ & V J 2 - ^atyJfy 

1 ~ 1 ~ l~ 

-^sd^psd^ps - -Stsd^ptsd' 1 pts - -o~vd I/ Jvd u J Vl , 

—5 t vd u Jf V d p 'JtVfi — -^tJv^ v {Jt VI x) — -jdtrJtv^ (JtTu^) 
+C P L c x (l) + C c T orr + C A , (14) 

r PC(l) _ fPC{l) fPC(l) XPCQ) . XPC0-) ( ,r\ 

L Lx ~ C S^S,Lx + C VW,Lx + C tS^tS,Lx + C tVHv,Lx { ib > 

1 -, v 
+ YQ a S\JTv f iJT + JtTufiJtr) 

+\[csd v {J Tv ,)d a {J7) + Ctsd^JtTv.mf^)], (16) 
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Table 3 

The different terms (except the scalar term) of the approximate Lagrange density. 

with primed coupling constants <Sj, 5i, 6{. These are not new parameters but 
can be related uniquely by a linear transformation to the given original couping 
parameters of the model 

{as,oits,oc v ,oitv,5s,5 t s,8v,8tv,0 T ,0 tv ) <= (a s ,a v ,atv,Ss,S v ,5 tv ) 

The actual transformation is a bit lengthy. It is given explicitely in appendix 
A.3. 

Finally, La is an electromagnetic Lagrangian density with an exchange cor- 
rection, it just remains to approximate the exchange part by the local density 
(Slater) approximation [30]. We have seen that in the Gordon decomposition 
representation, the dominant part of the £>j M and Cj densities can be replaced 
effectively by tensor terms and all possible mixings or combinations of non- 
linear terms in all orders. It means that if we perform Hartree calculations by 
taking into account tensor terms and all possible mixings or combinations of 
nonlinear terms of all orders, the effect of the dominant part of £> JM and Cj is 
already represented effectively, but still a part which manifests itself as an ef- 
fect beyond the Dirac equation and retardation effects are not yet represented. 
A strict Hartree theory including these additional effects can be constructed 
by starting off with the appropriate terms, but at the cost of introducing new 
parameters (the details about the mapping from Hartree-Fock to Hartree can 
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be seen in Appendix (3)). 



The new terms that emerge in our investigation lead to modifications of a) the 
spin-orbit force due to tensor terms, b) the isovector channel due to products 
of isoscalar and isovector terms and the inclusion of isovector- scalar terms 
and c) the density dependence of the model. They in turn lead to a different 
density dependence of the effective mass. A similar change of the spin-orbit 
force and the effective mass might allow a larger effective mass closer to the 
values predicted by Skyrme forces. It is interesting to discover the effect of the 
modified density dependence in both the isoscalar and the isovector channel. 
All these contributions arise without introducing new parameters. 

This unavoidably complicated mapping of the original exchange terms into a 
manageable Lagrangian in this rough approximation will now be applied to 
typical nuclear systems. 



5 Finite Nuclei 



Since the model discussed in this paper is meant to be applied to finite nuclei, 
observables connected with them should be considered when evaluating the 
consequences of our approximate treatment of exchange terms. We will take 
a look at a variety of different observables. The calculations are performed in 
spherical symmetry. 

As a preliminary step, we will apply further approximations. Let us take the 
following Lagrange density: 

r PC _ r , rPC , r PC(0) , r PC , p , r /-i 7 \ 

^NL off — ^-frec > *— l ^Lx cff *"NL ~r T i-Ax, < ) 

with 

£frec = E ~ m B )V a , 

a=A 

1 1 1 1 ~ 

~SvduJ^J Vli - l -8 tv dJ? v d v j tv ^ 
££?S = -\&xfts - ~ \o T J^{J Tvil ) - \e tT J^(J tTv ,), 

C A = -\d v A»d»A, + e E Va-fA^hl + r 3 )V a , 

1 a 1 

and 
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= —PspI - \lSP%. 

Thus now C^x^ ■> C^ rr , and C^Lx m Eq. (17) are neglected. CJ^ rr can be 
considered to give a small contribution, and it is assumed that parts of 
and ^nlx are already absorbed effectively in °f course, this assumption 
is not really true in the strict sense. We will study the following Lagrange 
densities for which we will determine parameter sets as described in the next 
section: 

(1) LH: Point-coupling model, no nonlinear terms, exchange omitted 

£LH = £frcc + £-1° + ^A- (18) 

(2) LHx : Point-coupling model, no nonlinear terms, approximate exchange 
included 

£lHx = £lH + ^Lx'eff + ^Ax- (19) 

(3) NLH : Point-coupling model, nonlinear terms, exchange omitted 

£nlh = £lh + C™. (20) 

(4) NLHx : Point-coupling model, nonlinear terms, nonlinear exchange omit- 
ted but linear exchange terms included 

£nlhx = £lh~x 

+ (21) 

Comparing these four versions should shed light on both the role of nonlinear 
terms and of exchange for realistic applications. 



5. 1 Determination of Parameters 



We follow the y 2 fitting procedure of Refs. [3,2], choosing the same set of phys- 
ical observables (binding energies, diffraction radii and surface thicknesses of 
16 0, 40 Ca, 48 Ca, 58 Ni, 90 Zr, lw Sn, 1M Sn, and 208 Pb ) and weights (0.2% rela- 
tive error for the binding energies, 0.5 % for the diffraction radii and 1.5% for 
the surface thicknesses) using the constant gap pairing correlation [3,2]. Our 
parameters are as, cty, a t y, 5s, $v and 5 t y (six parameters) for LHx and LH 
and additionally (3s and 75 for NLH and NLHx. 
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Table 4 

X 2 results from the models denned by Eqs. (18-21). Eb denotes the binding energies, 
Rd the diffraction radii and a the surface thicknesses 



Model 


a s 


ay 


a w 


Ss 


Sv 


$tv 


Ps 


75 


LHx 


27.632 


-14.881 




-0.517 


-0.152 


0.183 








LH 


-16.706 


12.829 


2.152 


-1.712 


0.670 


0.100 








NLHx 


-12.315 


2.193 


3.110 


-0.721 


-0.041 


0.045 


24.168 


-83.312 


NLH 


-14.726 


9.887 


1.345 


-0.746 


-0.058 


-0.377 


23.734 


-81.844 



Table 5 

Coupling constants of the various models 
5. 2 Results 



Table 4 collects the x 2 results from the four model variations. Here LHx, LH, 
NLHx, NLH denote parameter sets from Eqs. (18-21). E B , R d and a denote 
the binding energy, diffraction radius, and surface thickness. 

To see the difference among the four parameter sets clearly in binding energies 
of finite nuclei, we use the error in the binding energies (£e b ) an d the two- 
neutron (S2n) and the two-proton (£2?) separation energies from some isotopic 
and isotonic chains. The error in the binding energy is denned as 

c E B th — Eb ex P , 00 n 
c-Eb 7f ' \ ll > 

■CJB cxp 

while Si v and the are denned as 



S 2n = E B (N,Z)-E B (N-2,Z), 

S 2p = E B (N,Z)-E B (N,Z-2) (23) 

respectively, with Eb(N,Z) being the calculated binding energy. Fig. 1 shows 
that in the isotope chains, LHx has better £e b results than LH, but in the 
isotonic chains, LH performs better. Thus even though LHx has a better \\ 
than LH, it does not mean it must have better binding energy predictions for 
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Fig. 1. Error in the binding energies for selected isotopes (upper part) and isotones 
(lower part). 



all spherical nuclei. Although the Xe °f NLH and NLHx are not drastically 
different from that of LHx, they yield much better predictions than LHx for 
these isotopic and isotonic trends. We see here that the role of nonlinear terms 
for an acceptable £e b prediction is more important than that of the tensor 
terms, which in this case partly represents the exchange effect of the linear 
point-coupling model. NLH has a quality in Se b as good as NLHx. 

Fig. 2 shows that on the average, S 2n and S 2p from LHx are better than the 
ones from LH. Here we can see an important role of the nonlinear terms for the 
prediction of binding energies, because only NLH and NLHx can reproduce 
the S 2n and S 2p experimental data for almost all represented isotopes and 
isotones. The NLHx results are closer to the experimental data than NLH. 
The effect of the presence of tensor terms in the nonlinear model cannot be 
seen from £e b > but S 2n and S 2p show better trends. 
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Mass Number Mass Number 



Fig. 2. Two-neutron separation energies (S^n) for lead and tin isotopes (top) and 
two-proton separation energies (£271) for N=82 and N=126 isotone chains N=82 
(bottom). 

Errors in the diffraction radii (£R d )and the surface thicknesses (£ a ) are defined 
similarly as for the binding energies. The left picture of figure 3 shows that LH 
has an £ a of more than 20 % . LHx has an £ c value in the range below 20 %. It 
is still a bad result compared to NLHx and NLH which have an S a value below 
10 %. Although the xi °f NLHx is considerably better than for NLH, they yield 
similar predictions for the surface thicknesses. This confirms the result of Ref. 
[19], namely that the relativistic linear Hartree-Fock calculation still cannot 
give acceptable surface thickness predictions. Here it is clearly shown that the 
exchange effect has an important role in the surface thickness prediction. This 
figure also shows that only models which include the nonlinearities can give 
acceptable predictions. The right picture of figure 3 shows that the <?R d results 
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Fig. 3. Error in the surface thicknesses and diffraction radii for some selected iso- 
topes. 

of the parameter sets behave quite differently in each isotopic chain, but as the 
difference is not really significant, they may basically be considered to yield 
similar predictions in £R d - 

The definition of the error in spin-orbit splitting is the same as in the previ- 
ous cases. For the neutron case, all models give a bad prediction for 90 Zr(2p) 
and 48 Ca(ld). Only LH gives acceptable results in 208 Pb(3p), LHx and NLH 
in ie O(lp) and m Ca(ld), and NLHx in 208 Pb(3p), 16 0(lp) and m Ca{ld). For 
protons, LH cannot give acceptable results. LHx and NLH have three accept- 
able results while NLHx has four, all parameter sets cannot reproduce the 
experimental result of 48 Ca(ld) for protons. LHx is intermediate in quality 
between LH, on the one hand, and NLH and NLHx on the other for the 
single-particle spectra. From these facts it is clear that the exchange effects 
from the linear terms are important for the shell structure though quanti- 
tatively still insufficient, while the nonlinearity in the relativistic mean-field 
model clearly appears necessary for shell structure predictions. For the case 
of 208 Pb, LHx shows strong fluctuations in the surface part of the density and 
yields quite a large value of pp(r) in the center. This fact is probably due to the 
rough approximation used for the exchange of the derivative terms. We can 
see the role of the nonlinearity in NLHx which remedies the central density in 
comparison to LHx. 

The above results can be understood from the nonrelativistic analysis [32] 
as follows: the nonlinear models (NLH, NLHx) exhibit better predictions for 
finite nuclei observables than the linear ones (LH,LHx) because they have 
an adequat density dependence in the central and spin-orbit potentials. LHx 
shows better predictions than LHx in shell structure related observables be- 
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Fig. 4. Error in the spin-orbit splittings for some spherical nuclei (top) and the 
radial proton density and proton single-particle spectrum for 208 Pb (bottom). 

cause LHx contains a density dependence in the spin-orbit potential. LHx has 
better predictions than LH in binding energy and surface thickness, but the 
result is still far from acceptable (in contrast to NLH, NLHx) because LHx 
gets the density dependence in the central potential not from the dominant 
part C\ but from the minor part Cq. 

In summary, it is found that the exchange effect in the linear Hartree-Fock 
point-coupling model cannot be absorbed effectively into the coupling con- 
stants of a Hartree calculation. In agreement with Refs. [19,20] where the lin- 
ear Walecka model was used, a similar situation happens in the linear point- 
coupling model, namely that Hartree and Hartree-Fock calculations cannot 
give predictions which are close to experimental data. The nonlinear terms are 
needed to remedy this. Upon introducing the nonlinear terms, the exchange 
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effect from the linear terms does not drastically show up. Thus it seems that 
the nonlinear terms are more important than the exchange effect from the lin- 
ear terms. The validity of this presumption still needs to be checked, however, 
by the exact calculation, namely taking into account the exchange of both the 
linear and nonlinear terms exactly. 



6 Conclusion 

We have shown that formally, a Lagrangian densitiy C HF of the standard 
ansatz can be determined in a relativistic point-coupling Hartree-Fock sense 
that is not equivalent to that determined in a relativistic point-coupling strict 
Hartree sense, C H due to the exchange of the linear derivative terms and ex- 
change of the nonlinear terms. The equivalency can be obtained if we use a 
"complete" ansatz, but the cost we must pay is introducing more parameters. 
The exchange of the linear derivative terms is created by the densities Ci and 
£>j M . By using the Gordon decomposition, we can separate these densities into 
two parts. The role of the first part can be replaced effectively by tensor terms 
and all possible nonlinear terms, while the second part is a genuine feature of 
these densities. This second part is responsible for creating effects beyond the 
Dirac equation, nonlocal and retardation effects. An order of magnitude anal- 
ysis of the Ci and B,^ densities in the Gordon decomposition representation of 
these densities shows that these parts are expected to be small, yet the actual 
role of these parts in finite nuclei observables is still interesting to investigate. 
Also, as discussed before, taking into account the nonlinear exchange terms 
might alter the picture significantly. 

Our study, by using rough approximations in the Lagrangian densities, has 
shown gradual improvements in prediction for the finite nuclei. These results 
and the possibility that the effect of the exchange of the nonlinear terms 
and some neglected effects in the approximate forms of the densities Ci and 
Bin, which could reveal more improvements or even give a different picture 
in finite nuclei, support the importance of exact Hartree-Fock point-coupling 
calculations. 

Another interesting feature of such Hartree-Fock point-coupling calculations is 
the possibility to take into account terms which cannot appear in the Hartree 
calculation, for example the pionic degrees of freedom will contribute via the 
exchange terms. 

A reexamination of the ansatz is needed to improve the effect ivity of the model 
and a careful optimization procedure will be mandatory for fitting the model, 
both need elaborate and extensive work. For these reasons, we postpone the 
study of the complete quantitative Hartree-Fock model to future work. 
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A Appendix 



A.l From the Walecka to the Point- Coupling Model 

The advantages and disadvantages of both models in the context of the Hartree- 
Fock calculation are discussed. A Lagrangian density operator for the Walecka 
(W) model is chosen as in [21]: 

C W = CY + C% L (A.l) 

with 

+ E sAdJfd"^ - m^f^J - gMkM (A.2) 

i=S,V,R,A- 

and 

= - \hv'- (A.3) 



where the meson contents can be seen in table A.l. 



Meson 


K ( Field ) 


gi (Coupling Const) 


mi (Mass) 


(Coupling terms) 


Si 


a 


<p(x) 


9S 


m s 


1 


1 


id 




gv 


m v 




-1 


P 




9R 
2 


niR 




-1 


photon 




e 







-1 



Table A.l 



The meson contents of equation ( A.2). 

For the the fields used are ip(x), V^(x), R^x) and A^(x), denoting a 
scalar-isoscalar, a vector-isoscalar, a vector-isovector and the electromagnetic 
field operators, respectively. g t and mi are the coupling constants and masses 
of each field operator (f> ifl . We find that a difficulty appears to determine the 
scalar field operator due to the nonlinearities making it different from a simple 
Yukawa form. Some approximations have been done to overcome this prob- 
lem [22,23,24] but as far as we know there is not yet an exact Hartree-Fock 
calculation including mesonic nonlinear terms. Replacing 
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we can avoid the above problem and the standard definition for the linear field 
operators can be used as [21] 

Mx) = 4>^{x) +gi J d A yD(x - y, m^lr^, 

a 

*(x) =E(*«(a:)6 a + * a (x)dl), (A.4) 

a 

where ^ a (x), if? a and denote a nucleon, an antinucleon, and a meson 
wavefunction with ct enumerating the states. The operators a a j, b a and d a 
annihilate a free meson, a nucleon and an antinucleon with a momentum fc a , 
with the conjugate operators the corresponding creation operators. D(x,y) is 
the meson propagator which is defined as a solution of 



(d^ + m*)D(x - y, m 8 ) = 8\x - y). (A.5) 

The necessity to introduce an extra density dependence in Dirac-Brueckner- 
Hartree-Fock calculations that allows a simultaneous fit of the NN phase shift 
and the nuclear matter equilibrium point [25,26] provides a physical motiva- 
tion for this replacement. 

Further on, we can study the connection between the Walecka and point- 
coupling models by replacing the meson propagators with their low-order ex- 
pansion plus additional parameters ki as 



r d 4 n p-ip-( x ~y) i h. 

- / t^-A — rr — -v)- —^8\x - y ), (A.e) 

to arrive at the Lagrangian density operator of the point-coupling (RMF- 
PC) [17] model in section 1. 



A. 2 Gordon decomposition of the exact equation of motion 



After a straightforward Dirac matrix algebra calculation we have from the 
Gordon decomposition of the exact equation of motion that A 1 + A 2 + A 3 =0, 
with 
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+ \u;[{d^ a )^ a -^> a {d^ a )] 

+ l -W^d,[{d^ a )^ a - t a (3"tt a )], (A.7) 



A 2 = \d^ a a^ a ) - ^iU a ^(V a <7 a ^ a ) 

~\tW 2a d^d^ a a a ^ a ), (A.8) 



A 3 = m*y a ^ a + V^ a V a . (A.9) 



Let us define m* = m* + Am* and V a = V a + AV a , where Am* and AT/ Q are 
the parts of m* and V a of order t> 2 . Note that up to this point the treatment of 
exchange is still exact. Now, if we apply the order-of-magnitude estimation of 
section 3 to study the magnitude in every term in A\ through A 3 and assume 
that the terms of order > v 2 are small and can be neglected, it is obvious 
that only the first term in Ai and A 2 will survive and the Am* and AV a 
contributions vanish in A 3 . Therefore in this limit we will obtain Eq. (6) of 
section (3). 
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A. 3 The coupling constants of the approximate Lagrangian 



1 4 12 1 4 12 

a s = a s - (-a s + -a v + —a w ) - 2m B (-5 s + -S v + yW)], (A. 10a) 

t 1 2 6 , o2/^r 2 r ^ r \ 

= «v + (-g«5 + gtty + g«ty) + 2m B (--^ 5 + -d v + -6 tv ) 

+2m%(^6 s + p v + ^-S tv ) (A.lOb) 

14 4 14 4 

cits = (-g«s - gtty + gtttv) + 2m B (--5 s - -5 V + -5 tv )] (A.lOc) 

12 2 1 2 2 

a t v = a t v + (- g«5 + g"v - g«ty) + 2m B (--5 s + -<V - -<W) 

-2m%(~S s - p v + ptv)] (A.lOd) 
1 4 12 

5 S = 5 S + (y^s + — 5 V + — <V) (A.lOe) 

Sv = Sv - (--U? + ^-Jv + ^<W) (A.lOf) 
lb lb 32 

fits = + ^6 V - -^8 tv ) (A.lOg) 

<V = <V ~ (~Y^ S s + ^S v - ^5 tv ) (A.lOh) 

9 T = Am B {-^5 s + ^5 V + ^<V) (A.lOi) 

9 tT = 4m B (^-5 s + ^-5 V - ^-5 tv ) (A.lOj) 
lb lb lb 



.1 4 . 12 r ^ 

cs = ( r6 5s + -5 v + -6tv), 

c v = (-Y6 5s + Y6 5v+ 32 5tv) ' 
cts = (-^5s + ^5v-^5tv), 

ctv^i-^s-^v + ^tv), 



(A.lOk) 



A complete mapping from Hartree-Fock to Hartree 



It has been shown that the usual ansatz for the point-coupling RMF model, 
when considered within the Hartree-Fock approximation, leads to equations 
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that deviate from the structure of a Dirac-equation. The source for these 
differences are the exchange terms with derivatives, leading to densities that 
do not occur in the standard Dirac-Hartree model. 



It should be mentioned, however, that the reason for these deviations is the 
form of the Lagrange density operator, with which we start off. We could as 
well choose a Lagrangian meant for a strict Hartree treatment, which would 
lead to the same structure as we obtained it in this investigation. A simi- 
lar situation holds for the Skyrme-Hartree-Fock theory. The mapping can be 
achieved by adding terms that lead to the new densities in the Hartree-Fock 
case. The Lagrangian density would have to be complemented by the following 
terms: 



derivative _ 



(A.11) 



here etj, f3 3 and 9 3 are free parameters. This would include, however, the intro- 
duction of additional parameters, whereas in the HF case, no new parameters 
arise. The ansatz complemented by Eq. (A. 11) leads to the same type of equa- 
tion in both the Hartree and the Hartree-Fock approximations. 

Though the mapping of the Hartree-Fock (and the corresponding effective 
Hartree) theory to a related strict Hartree theory can be done in principle, it 
would be at the cost of introducing and motivating new terms and additional 
coupling constants. If we also consider higher order terms [27], a complete 
mapping to a Hartree theory would introduce many new parameters. Experi- 
ence will tell which ansatz for the effective interaction will prove to be most 
adequate for a Hartree-Fock theory for finite nuclei. 
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